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We discuss the local density approximation approach to calculating the ground state energy 
of a one-dimensional Fermi gas containing a single impurity, and compare the results with exact 
numerical values that we have for up to 11 particles for general interaction strengths and up to 30 
particles in the strongly interacting case. We also calculate the contact coefficient in the strongly 
interacting regime. The different theoretical predictions are compared to recent experimental results 
with few-atom systems. Firstly, we find that the local density approximation suffers from great 
ambiguity in the few-atom regime, yet it works surprisingly well for some models. Secondly, we find 
that the strong interaction theories quickly break down when the number of particles increase or 
the interaction strength decreases. 
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I. INTRODUCTION 


As cold atomic gas capabilities advance Hi, a large 
number of interesting quantum systems can now be built 
using the cold atom toolbox iJi. The structure and 
dynamics of one-dimensional quantum systems is a spe¬ 
cific venue in which cold atoms have driven a new surge 
of excitement as a number of important models of many- 
body quantum physics can now be probed experimen¬ 
tally [6l-[ll|. This includes the famous strongly interact- 
ing Bose systems studied by Tonks [l^ and Girardeau 
[iq . More recently, it has become possible to push these 
studies to smaller particle numbers so that few-body 
physics in one-dimensional traps with controllable inter¬ 
actions can now be realized. A number of such studies 
in recent years using two-component Fermi systems have 
probed such concepts as fermionization [l^ , pairing d, 
polarons , and the Hubbard [2lj and Heisenberg mod¬ 
els [12, all in the limit where the system contains just a 
few particles. These developments have sparked a lot of 
interest in the few-body community and a number of new 
theoretical studies looking at various aspects of these two- 
component Fermi systems in one dimension have emerged 

[Hil. 

In the present paper we are interested in the recent ex¬ 
periments that have shown how one may build a Fermi 
system consisting of an impurity and a Fermi sea one par¬ 
ticle at a time [12 . This amounts to building a so-called 
Fermi polaron from a few-body perspective. The experi¬ 
mental results presented in Ref. [12 suggest that already 
for a system consisting of only six particles (five identi¬ 
cal fermions and a single impurity of the same mass as 
the other particles) one can see the emergence of a many- 
body system taking place. This is a rather remarkable 
result and in the present paper we would like to explore 
this in more detail. The experimental results consists 
of a series of measurements of the energy for different 
interaction strengths and for different system sizes and 
have previously been compared to theoretical predictions 
in the literature [12. However, there are some ambigu¬ 


ities in this comparison arising from the fact that the 
experimental system is in a trap while a lot of theoreti¬ 
cal work considers a homogeneous system. Proposals for 
merging the homogeneous results with a trapped system 
are known and one of our goals here is to assess how well 
they work as a function of particle number and interac¬ 
tion strength. We will take advantage of the fact that we 
have access to both very accurate numerical calculations 
for all interaction strengths for up to eleven particles and 
in the strongly interacting regime for up to thirty par¬ 
ticles in a harmonic trap. Using these results we will 
shed some light on the questions concerning how many 
particles constitutes ‘many’ in one-dimensional impurity 
systems and when such systems can be considered in the 
strongly interacting regime. 


II. MODEL 


We consider a system of N spin-i fermions with mass 
m confined in a one-dimensional (ID) harmonic trap with 
oscillator frequency w. The fermions in the gas are spin- 
polarized with N — \ particles in, say, the spin-up state, 
while the remaining particle (the impurity) is in the spin- 
down state. The particles interact via a short-range inter¬ 
action modelled by a delta function of coupling strength 
g. The system is described by the Hamiltonian 


N 


^ = E 


1 2 2 


Pg'^Sixi-Xj) . ( 1 ) 


Kj 


In the interaction term, the sum runs over all pairs of 
particles, but due to the antisymmetry of the total wave 
function under exchange of any two majority particles, 
the majority particles do not couple to each other. Thus 
the interaction term above is equivalent to only summing 
over the — 1 pairs consisting of the impurity and each 
of the majority particles. As has been recently demon¬ 
strated experimentally, the interaction strength may be 
varied from small to large positive and negative values 
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We will limit ourselves to a repulsive interaction, 
g > 0. Throughout this paper, we use harmonic oscilla¬ 
tor units where oj = m = h = 1. 


III. STRONG COUPLING REGIME 


We start by discussing the Tonks-Girardeau (TG) limit 
in which the interaction between the particles is infinitely 
strong, 1/g —>■ 0. In this limit, the ground state becomes 
N\/{{N — 1)! • I!) = Wfold degenerate with the energy 
= iV^/2, and the wave function vanish whenever any 
two particles meet. The superscript t denotes the fact 
that the system is trapped, and the subscript oo refers to 
the interaction being infinitely strong. Notice that the 
degenerate energy in the ground state manifold is also 
the ground state energy of an ideal non-interacting Fermi 
gas of N particles in a harmonic trap, i.e. the particles 
occupy the N lowest-lying energy states in the trap. 

The degeneracy is lifted by moving slightly away from 
1/g = 0, and the wave function is lifted slightly away 
from zero at the surfaces where the impurity meets a 
majority particle. This allows the impurity to switch 
position with its neighbors, and we can therefore think of 
the particles as occupying sites on a lattice. The system 
can be described by an effective Heisenberg spin chain 
model as discussed in [H, [13, IH, IHl • To linear order in 
1 / (/ <C I, we may cast the Hamiltonian as 


i? = i ^ ^ [I - . ^^-+1] , (2) 

^ fc=l 

where cr^ = (ct^, ctJ) are the Pauli matrices acting on 
the spin of the particle at site k, and ak is a geometric 
coefficient determined by the trap potential. These are 
the same a-coefhcients discussed in [i^, where the case 
of a general confining potential, V(x)^ and an arbitrary 
number of spin-up and spin-down particles is considered. 
The general case captures the behavior of both bosonic 
and fermionic particles near the TG limit. In Ref. [56l| it 
is shown that the geometric coefficient ak generally can 
be expressed as: 


N N N-l-k 

i=l j=l 1=0 


(^-ly+j+N-k _ I _2 

k-1 


n 




X J dx‘^{V{x) - 


^det [(i?(x)-AI)(^^)] 
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dx 


( 3 ) 


- x—b 


To compute this we need access to the N lowest- 
energy single-particle wavefunctions 'ipi{x) solving the 
Schrbdinger equation + V (a:)] = Eiijji. In the 


present case, this is the analytically solved Schrddinger 
equation for a single particle in a harmonic oscillator. In 
Eq. ([3]) i?(x) is a iV X iV symmetric matrix with the 
mn’th entry given by dyi>m{y)4’n{y), and ( de¬ 
fines a minor obtained by removing the I’th column and 
the j’th row. For most potentials such as smooth poten¬ 
tials like the present case, the integration limits will be 
6 = 00 = —a, and for some non-analytic potentials such 
as a hard box, the integration will be confined to some fi¬ 
nite region of space. For symmetric confining potentials, 
as the one studied in this paper, the coefficients are also 
symmetric, ak = a^-k- 

The main motivation for the present work is that we 
are able to determine these geometric a-coefficients ex¬ 
actly in the general case of an arbitrary potential, see 
Ref. [ 5 ^ for a detailed presentation of our method in¬ 
cluding a derivation of Eq. Q. The method is imple¬ 
mented as a highly efficient numerical program released 
as open source software [s^. We stress the importance 
of our capability to perform these calculations numeri¬ 
cally exact; we do not use Monte Carlo integration or 
any other approximate method. For this paper, we cal¬ 
culate the geometric coefficients numerically exact for a 
harmonic potential for TV < 30. The coefficients are given 
in Section [X] Calculating the coefficients for such a high 
number of particles as 30 is an extremely demanding com¬ 
putational task. However, the method we have developed 
enables us to perform these calculations on a small lap¬ 
top to a very high level of precision in a matter of hours 
for 30 particles. For a few particles, the calculations are 
done almost immediately. 

Since the harmonic potential is such a well-studied con¬ 
fining potential, this allows us to compare our exact re¬ 
sults with the approximate methods studied in the lit¬ 
erature. In particular, we are interested in comparing 
methods for calculating the ground state energy in the 
vicinity of the TG limit. To linear order in 1/g, we write 
the ground state energy as 

El = El,--C, (4) 

9 

where C is the eigenvalue of the sum in Eq. © for the 
ground state, i.e. the largest eigenvalue. We denote C the 
contact coefficient [^, . Access to the exact geometric 

coefficients ai allows us to diagonalize the Hamiltonian in 
Eq. (I2|) and calculate the contact coefficient exactly. For 
TV < 30 we compute C and want to compare our results 
with previous approximate formulas for this constant. 


A. Previous results 

Now we will review some previous calculations of the 
contact coefficient C in Eq. 0. However, since the pre¬ 
vious theories may not take the trap into account, we 
will have to adapt the existing formulas in order to di¬ 
rectly compare them with our results. In the following. 
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we will distinguish between the free system and the (har¬ 
monically) trapped system. Furthermore, we will discuss 
different regimes of the interaction strength g and the 
particle number N. 

A pioneering piece of work was done by McGuire [52| 
who calculated the ground state energy for a free ID 
Fermi gas with a single impurity for an arbitrary in¬ 
teraction strength g subject to periodic boundary con¬ 
ditions. Thus the trap needs to be taken into account. 
This problem was recently addressed by Astrakharchik 
and Brouzos who applied the local density approxima¬ 
tion (LDA) to McGuire’s result in order to obtain an 
expression for a harmonically trapped system [^ . How¬ 
ever, McGuire assumes that the number of particles N 
is even and large. This subtlety does not appear to be 
explicitly discussed in Ref. [s^. Namely, the implicit con¬ 
dition that N is very large makes it interesting to com¬ 
pare the result with our exact results for small N. To 
get the desired contact coefficient C, we should expand 
Astrakharchik and Brouzos’ LDA result around 1/g = 0 
to first order and read off the coefficient. 

McGuire’s energy expression was found for a free inter¬ 
acting system subject to periodic boundary conditions 
(PBG). It would be interesting to consider a hard wall 
boundary condition (HWG), too. This problem was 
solved by Oelkers et al. for a free gas in the strong cou¬ 
pling regime which is exactly the regime we are inter¬ 
ested in. In fact Oelkers et al. found results for both PBG 
and HWG for finite (even) N, and not just in the limit of 
large N. Just as Astrakharchik and Brouzos applied LDA 
to McGuire’s free result, we shall apply LDA to Oelkers’ 
free results. Since the LDA method maps the free sys¬ 
tem to a trapped system, we would intuitively think that 
HWG would produce a better final result, since a hard 
wall would resemble a harmonic potential more that a 
periodic boundary. This will turn out not to be the case, 
but as we shall also see, these LDA results are subject to 
high uncertainties in the small N regime. Nevertheless 
this is an important regime as current experiments study 
the behaviour of few particles. 

The starting point for our analysis are the results sum¬ 
marized in Table [J Here we have recorded whether the 
energy is for a harmonically trapped or a free system, 
which boundary condition is used and which values of 
the interaction strength and particle numbers it is valid 
for. 


Reference 

Trapped 

BC 

5 

N 

McGuire [M] 

No 

BBC 

any 

A —>■ oo 

Astra. [30] 

Yes 

PBG 

any 

A —>■ oo 

Oelkers [^] 

No 

PBG 

5 -)■ oo 

even 

Oelkers 

No 

HWG 

g ^ oo 

even 


TABLE I: Summary of some previous results. 


E^f = Wf? 



5 = 0 1/5 ~ 0 1/5 = 0 


FIG. 1: Illustration of the local density approximation. The 
top part of the picture shows the free system solved for some 
boundary condition (green), whereas the bottom part shows 
the trapped system (red). As a way of representing an inter¬ 
acting sea of particles, we have drawn waves on the surface 
with increasing sharpness for increasing interaction strength. 
This should only be understood as an intuitive guide to the 
eye. The Fermi level is defined in the non-interacting case and 
denoted by an arrow, while the total energy of the system is 
denoted by a circle. 


B. The local density approximation 

The strategy to map a result from the free case to 
the trapped case is by applying LDA. Therefore, we will 
briefly discuss LDA, and most importantly we will en¬ 
counter a weakness in the method for small values of N. 

Denote by (Eq) the energy of the free (trapped), 
non-interacting system. Suppose that we know the en¬ 
ergy states of these two systems, and that we therefore 
can associate with them Fermi levels and Fermi momenta, 
denoted Ep (Ep) and kp (kp), respectively. Suppose fur¬ 
thermore that we know an expression for the energy of 
the interacting system, but only in the free case. Denote 
this energy by 


El = Ei + /^El(ky), (5) 

where is the correction due to interaction, which 

will generally depend on kfp. Now, the quantity sought 
after is the energy for the interacting and trapped system. 
Within LDA, this is found by mapping 

Eq Eq and kpi-^kp. (6) 

The interaction correction in the trapped case retain its 
functional form from the free case, but is now evaluated 
at the trapped Fermi momentum instead of the free Fermi 
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FIG. 2: The Fermi level, should lie between the highest energy of the occupied states and the lowest energy of the 
unoccupied states for the non-interacting system. In this limit the N — 1 states with lowest energy are filled up. Notice that 
in PBC for N odd, the highest occupied state and the lowest unoccupied state share the same energy, and thus the Fermi level 
is well-defined, but that we have a certain freedom in the other cases. This freedoms leads us to define the tuning parameters 
APbc,aHwc^;,ho^ 


momentum. Thus the LDA expression for the interacting 
and trapped system is 

El^El + ^El{k^p) . (7) 

An illustrative sketch of the LDA procedure is shown on 
Figure [T] 

We have just described the LDA method if our starting 
point was the known non-interacting limit for both the 
free and trapped system. If, however, we approached a 
finite interaction starting from the infinitely strong inter¬ 
acting case, l/g = 0, we would instead have mapped 

^oo kpi-^kp, ( 8 ) 

where {E^^) is the energy of the free (trapped) in¬ 
finitely strongly interacting system, both assumed to be 
known. We use this approach in Subsection IIII PI 

The employment of LDA relies on the knowledge of 
the Fermi momenta, k^p, for the free and trapped (non¬ 
interacting) systems, or equivalently the Fermi levels, 
Ep. Now the problem is that the Fermi level (also called 
the chemical potential) is not well-defined if the systems 
contains a finite number of particles. We know that the 
Fermi level should lie between the highest energy of the 
occupied states and the lowest energy of the unoccupied 
states. But exactly where between these two energies is 
irrelevant for the occupancy of the states. The problem 
is illustrated on Figure [5] for the four relevant situations. 
In the non-interacting limit, the N — 1 identical fermions 
hll up the — 1 states of lowest energy. Recall that 
the momentum for a particle subject to PBC, HWC or 
in a harmonic oscillator is given as 


,PBC 27r 

" ~ L ^ 

* = 0, il, i2. 

,HWC _ 

f = 1,2,3,... 

* L ’ 

=6-1^2^, 

* = 0,1,2,... 


where L is the length of the system and b = y'h/rriLo is 
the standard oscillator length, 6 = 1 in our units. Then 


the energy at the Fermi levels can be described using the 
following expressions for the Fermi momenta: 

f (Af - 2 -f if N is even 

fiN-l) if Af is odd ’ 

fcHWC = 1 + yHWC) ^ 

Lj 

kf^ = 6-y2(N- I -b AHO) . 

Here A“ G [0,1] with a = PBC, HWC, HO is some tuning 
parameter that allows us to probe the energies between 
that of the highest occupied state and the lowest unoc¬ 
cupied state. Picking A“ = 0 corresponds to picking the 
Fermi level at the highest occupied state and A“ = 1 
corresponds to taking the Fermi level at the lowest un¬ 
occupied state. Another appealing choice is A“ = 1/2, 
corresponding to taking the Fermi level right in the mid¬ 
dle. For the time being, we will carry around the tuning 
parameters, but at some point we would like to pick spe¬ 
cific values. 

Note in particular that the ambiguity of the Fermi mo¬ 
mentum disappears in the thermodynamic limit, because 
the difference between the highest occupied state and low¬ 
est unoccupied state becomes insignificant when N ^ oo. 
This is reflected in the fact that the arnbiguity of the 
Fermi momentum is not discussed in Ref. [30| , where also 
the concepts of the Fermi level and Fermi energy, which 
only equal each other in the thermodynamic limit, are 
used somewhat interchangeably. But if we want to de¬ 
rive energy expressions that apply to finite N, we should 
be careful when choosing the Fermi momentum. 


C. McGuire, Astrakharchik and Brouzos (PBC, 
thermodynamic limit) 

In Ref. McGuire finds the following expression for 
the ground state energy of the interacting, free system 
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subject to PBC: 


AEl = 


f _ 7 


fl-7 + 

A + 

27^^ 

tan-3 ^ 

4 

^27r 

7 y 

27r 


(9) 

where 7 = gn/kp^^. The result applies to all values 
of the interaction strength g, but in deriving Eq. Q 
McGuire converts a sum to an integral letting N ^ 00 
and L ^ 00 with the density N/L held constant. This 
sum to integral conversion can be done in several ways 
introducing some degree of freedom in consistent 

with the discussion in the previous section. Notice that 
picking = 1/2 would set k^^^ = {N— 1)^1 L for all 
N. This would imply a vanishing interaction correction 
for 1 as it should. [ 6 l| 

We now sketch how Astrakharchik and Brouzos imple¬ 
ment LDA on McGuire’s free energy exp ression @ to 
find the energy of the trapped system The LDA 

expression for the energy of the trapped system © yields 


El ^ El, + AAf (fcHO) 

(iV- 1)2 -h 1 
“ 2 
(A:HO)2y 
2 ^2 

( 10 ) 


1-7 + 

^7‘ , 

27r/ 

tan-3 TL 

4 

i 27r 

7* , 

27r 


with 7 * = gn/k^‘^. As before, we have to pick = 
1/2 and thus = \/2{N — 1) in order to ensure that 
the energy correction vanishes for A —>■ 1. Since we are 
interested in the energy in the case of strong interaction, 
we expand the above general expression to first order in 

1 / 5 : 


iV2 1 

“2 


where the desired contact coefficient is given as 


( 11 ) 


C =—(iV-1)3/2 . (^2) 

Stt 

The above LDA expression is directly comparable to our 
exact results. Notice that the LDA method introduced 
ambiguities in the choice of the Fermi momenta, but that 
known physics could be used to restrict the choice and get 
an unambiguous final result. Unfortunately, this cannot 
be done in the calculations based on Oelkers’ free results 
for a finite N in PBG and HWC. 


D. Oelkers (PBC and HWC, finite particle 
numbers) 

Before we compare the approximated expression in 
Eq. (fT^ with our exact calculations, we wish to derive 
two other approximated expression for the contact coef¬ 
ficient. These two expressions will rely on the results de¬ 
rived by Oelkers et al. in Ref. . In this reference, the 


authors calculate the energy of a free strongly interact¬ 
ing system using PBC and HWC. Contrary to McGuire 
who assumed large iV, these results are valid for finite 
(but even) N, which is what we ultimately are interested 
in. Our plan is now to apply the LDA method in order 
to get expressions for the trapped systems. We start by 
considering the case of PBC. 

Periodic boundary condition. To first order in 1/g, we 
find the free ground state energy to bej^ 


E^ = - (N^ 


■2N) 


= Ei 




^ gL. 
■2N). 


(13) 


We now follow the mapping prescription ([5]) and map 
EL ^ = /V2/2 and k^^ ^ feg®. Using the ex¬ 

pressions for k^^^ and k^^, we can rewrite the mapping 
as 


L 7r(A - 2 + 2APBC) ’ CU 

where we have kept the tuning parameters unspecified. 
Applying the above substitution to the result in (fT^ , we 
obtain the energy of the trapped system within LDA: 


El 


iV2 


1 8^2 

g Stt 


(A3 + 2N) 


(iV - I + AH0)3/2 
{N-2 + 2APPC)3 


(15) 


In the thermodynamic limit, N ^ 00 , the contact coef¬ 
ficient read off Eq. (1151) and the one in Eq. (fT^ tend to¬ 
wards the same asymptote given by Coo = 8-\/2 /(Stt)^ 3/2, 
regardless of the choice of the tuning parameters. This 
is indeed reassuring, but we are mostly interested in re¬ 
sults for small N. Eq. (fT5l) is derived for even N, but let 
us extend the domain also to odd N in order to restrict 
the tuning parameters by requiring the vanishing of the 
contact coefficient for A = 1. Since we also require C to 
be non-divergent for finite A, we end up with A^® = 1/2 
and APP® € ]0; l/2[U]l/2; 1]. While this fixes A^®, there 
is still ambiguity in AP®®. Choosing AP®® = 1 yields 


C = 


8^/23 , oAn(^- 1)^^^ 


Stt 


+ 2N)- 


N3 


(16) 


Hard wall condition. The ground state energy with 
HWC is found in Ref. 


Ei = 


12L2 


(2^3 - 3^2 -h N) 


1 8 2 / ^ A 

1 -r cos 

gL \2N J 


(rq 

This result is consistent with the energy found in Ref. 
in the thermodynamic limit {N,L —>■ 00 and N/L con¬ 
stant). In a completely similar way as before, mapping 
!->• / 2 and >->• k^^ yields the energy 
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FIG. 3: Comparing the approximated contact coefficients 
based on the results of Oelkers et al. [s^ and McGuire [5^ 
with our exact results, scaled with the common asymptote for 
N ^ 00 . For Oelkers HWC we also change the tuning param¬ 
eters to see how that influences the result. The dotted lines 
show the choices A**° = 1 and = 1/4 (upper curve) and 

A^° = 1/2 and = 1 (lower curve), where we for clarity 

have taken to be a continuous parameter. 


of the trapped system: 


=— - -^(N^ -^N^ + iiV) 

2 g Sir ^ ^ ^ ^ ^ 

(A1-| + AH0)3/2 / tt \ 

"" (iV - 1 + AHWC)3 cos ( 2 ^) . (18) 

As expected, the contact coefficient goes towards the com¬ 
mon asymptote Coo in the limit N ^ 00 . Again we wish 
to extend the result to odd N, not just even N. To en¬ 
sure that the contact coefficient is real and well-defined 
for any integer A^ > 1 and zero at A^ = 1, we should pick 
either A^® = 1/2 and A^'^ci = 0 or A^c> g [1/2; 1] and 
AHWC g]0;l]. For the choices A^^* = _ \j2^ we 

get 


C = 




|Af2 + lAf) 


(Ai- 1)3/2 


cos^ 



In deriving Eqs. (ffSl) and di), we had to specify values 
for the tuning parameters, but notice that we could just 
as well have chosen other values within certain bounds. 


E. Comparing contact coefficients 

We now wish to compare the contact coefficients that 
we derived in the previous two subsections. First of all, 


we want to see, whether the approximated results com¬ 
pare well with our exact results. Secondly, we want to 
investigate how the choice of the tuning parameters inter¬ 
feres with the result for small values of N. On Figure [3] 
we have shown our exact results compared with the three 
LDA expressions in Eqs. (fT3]| . (fTHll and (fT31) . 

We see that the LDA results tend to undershoot com¬ 
pared to our exact results. We also wanted to compare 
HWC results and PBC results, because our intuition told 
us that a hard wall box ‘looked more’ like a harmonic trap 
than a free periodic potential did. Contrary to our ex¬ 
pectation, the HWC result does not appear to be better 
than the PBC results. On the other hand, the ambiguity 
in the LDA results introduced by the choice of the Fermi 
level and manifested in the tuning parameters makes it 
impossible to compare HWC and PBC in a unique way: 
We have to choose values for the tuning parameters. On 
Figure [3] we also show the HWC result for other values 
of the tuning parameters, resulting in drastically altered 
contact coefficients. Even for TV ~ 30, the choice of the 
tuning parameters matters, but eventually it will become 
insignificant as TV becomes larger and larger. 

In our opinion this constitute the most important les¬ 
son about the LDA method: For values of TV that realisti¬ 
cally can be probed experimentally today or in the com¬ 
ing years, the partially arbitrary position of the Fermi 
level has a very large effect on the LDA predicted con¬ 
tact coefficient. In the light of this conclusion, it seems 
like an improbable stroke of luck that McGuire’s result 
found in the limit TV —> 00 works so well when extrapo¬ 
lated to the finite TV regime. 


IV. COMPARING MODELS TO EXPERIMENT 

The ground state energy for the impurity system was 
found experimentally in Ref. for TV = 2,... ,6 and 
g = 0.36, 1.14 and 2.80. Naturally, we want to examine 
how the theories that assume large TV and/or g compare 
with each other and the experimental results. 

McGuire assumes that TV —>■ 00 , so we will use Eqs. (iini) 
and m as the approximated theoretical energies assum¬ 
ing that g is finite and very large, respectively. We use 
our exact contact coefficients as a theoretical prediction 
assuming finite TV, but with g being very large. We also 
have exact numerical calculations of the energy for finite 
TV and finite g, but this is a very demanding computa¬ 
tion, so we only have reliable results for TV < 11 (using 
the effective interaction method recently introduced to 
address Fermi and Bose systems [ 53 ). These four 
theoretical predictions exhaust the four combinations of 
finite/infinite g and TV. We compare them with the ex¬ 
perimental data on on Figured 

Not surprisingly, the most general theory obtained 
from full numerical calculations for finite values of g and 
TV is in very good agreement with the experimental data. 
But we also see that the theory assuming finite g and 
TV —> 00 works equally well for very small values of TV. 
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FIG. 4: Interaction correction to the energy of the non-interacting system as function of the total number of particles. The 
colors refer to different values of the interaction strength. On the top figure, we compare experimental data with theories for 
finite g. For < 11 we have reliable results from a full numerical treatment for finite g and finite N (+). Notice that the 
agreement with the McGuire formula for finite g and oo in Eq. m (—) is excellent in the few-body regime and is 

expected to improve for increasing N. On the bottom figure, we compare experimental data with theories assuming g —>■ oo. 
We compare the results for finite N using Eq. and the exact geometric coefficients in Appendix 0 ( • ) and the McGuire 
formula in the limit —>■ oo in Eq. m (•••)• To illustrate how the g ^ oo theories fall off for increasing N we also show a 
piece of the finite g McGuire formula for 5 = 10 on the bottom plot. 


This is indeed an unexpected result, especially when we 
recall from the previous section that the exact position 
of the Fermi level chosen in the LDA could have a huge 
impact on the final result for small N. 

Let us now examine the theories assuming g ^ 00 . 
We immediately see that our exact result and the LDA 
result are in very good agreement as a direct consequence 
of the agreement between the exact contact coefficients 
and the McGuire result in Figure[3l We know that the in¬ 
teraction correction goes like Ai?* ^ —N^^'^/g for large 
N, so for suitable large N and/or small g the theory 
must break down (at some point the interaction correc¬ 
tion even becomes negative). This is clearly seen when 
comparing the top and bottom part of Figure [D It is 
seen that strongly interacting theories deviate from the 
general theories when N increases or g decreases. The 
interesting question is then for which N and g does the 
strongly interaction theory describe nature? The largest 


experimental interaction strength probed in Ref. is 
g = 2.80, but as is seen on the bottom part of Figure® 
theory and experiment only matches for N = 2, so this 
is clearly not very strong interaction. For g = 10, we get 
a good agreement between the theories assuming finite g 
and large g for N < 12, but then the theories split up. 
We find it surprising how quickly the strong interaction 
theory breaks down. 

V. SUMMARY 

We have considered a trapped one-dimensional Fermi 
system interacting with a single impurity of the same 
mass from the point of view of energetics. We have dis¬ 
cussed how the local density approximation served as a 
way to map the ground state energy of the free system 
into the ground state energy of the trapped system. We 











also saw that the method could introduce ambiguous re¬ 
sults for finite N that heavily influenced the predicted 
energies. This naturally lead to a skepticism for LDA 
results for small N, which are currently used to describe 
experiments [ 2 ^. Despite the erratic nature of the LDA 
results, the energy expression based on McGuire’s result 
works surprisingly wellj^. [^. 

We also discussed how the strong interaction limit must 
be approached with caution as the particle number in¬ 
crease. Comparing theories for finite interaction and 
strong interaction, we found that g = 2.80 cannot be 
classified as strong interaction, and for g = 10 the theo¬ 
ries start deviating significantly already around A^ ~ 12. 
This is consistent with the observation that the natural 
interaction strength in the harmonic trap is proportional 
to g/'/N [ 3 ^ (see also the earlier discussion in Ref. [13 )■ 

Remark. During preparation of the revised version of 
this manuscript, we became aware of a paper by Deuret- 
zbacher et al. which describes a method for com¬ 
puting the geometric coefficients and presents results for 
harmonic traps with up to 30 particles using a method 
that has some common features with the computational 
approach used here. 
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Appendix A: Geometric coefficients for a harmonic 

trap 

Here we tabulate the geometric coefficients, Oi, in 
Eq. ([2|) in harmonic oscillator units. For completeness we 
note that the explicit unit is {hu})‘^b, where b = \Jh/mu} 
is the standard oscillator length. To the right of every 
indicated value of N we record ai, 02 , • ■ • up to ajsi /2 
for N even and q;(jv-i )/2 for N odd, continuing on the 
next line if necessary. The remaining geometric coeffi¬ 
cients are given by the symmetry property = a^-k 
for k = 1,...,A — las there are N — 1 coefficients for each 
N. We have calculated all A — 1 coefficients for each N 
and used the symmetry property to estimate the numeri¬ 
cal precision on the coefficients by noting the digits in au 
in agreement with the digits in a^-k- In this appendix 
we have at most given 10 significant figures even though 
many of our calculations are even more precise. Gen¬ 
erally a higher precision can be achieved for the price of 
higher calculation time. The coefficients given here agree 
with previous calculations which have been done for 10 
particles or less [H, 113,1131 . 
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N ai = ajv-i 02 = OAr -2 

2 \/^ = 

0.7978845608 

3 3V(250F) = 

1.346430196 

4 1.787645708 2.346508058 

5 2.166057718 3.177197531 


6 

2.5021784 

3.902098540 

4.357116131 






7 

2.807397825 

4.552904442 

5.400326410 






8 

3.088795182 

5.148133894 

6.345017625 

6.738585693 





9 

3.351118130 

5.69971481 

7.214307411 

7.959044430 





10 

3.597730090 

6.215859197 

8.023484036 

9.08789432 

9.439679835 




11 

3.831114914 

6.702502948 

8.78328981 

10.14280881 

10.81522991 




12 

4.053167612 

7.164099681 

9.501616401 

11.13640807 

12.10509843 

12.42610595 



13 

4.265372369 

7.604090937 

10.18446319 

12.078086 

13.32332240 

13.94120413 



14 

4.468916965 

8.025200651 

10.83651694 

12.97506385 

14.48046497 

15.37564625 

15.67276485 


15 

4.664769312 

8.429628019 

11.46152119 

13.83303366 

15.584720 

16.74091572 

17.31574325 


16 

4.853730361 

8.819178529 

12.06252191 

14.65657900 

16.64259962 

18.04599032 

18.88235477 

19.16023435 

17 

5.036471720 

9.195355790 

12.64203656 

15.44945457 

17.65938206 

19.2980535 

20.38221388 

20.92192258 

18 

5.213563048 

9.559427571 

13.20217385 

16.21478205 

18.63941758 

20.50296313 

21.82305903 

22.61094361 


22.87289739 








19 

5.385492433 

9.912474388 

13.74472084 

16.95518993 

19.58634200 

21.66557213 

23.211233 

24.23549760 


24.74584444 








20 

5.552681814 

10.25542595 

14.27120748 

17.67291605 

20.50323106 

22.78995562 

24.55201753 

25.80233751 


26.54935034 

26.79783773 







21 

5.715498860 

10.58908900 

14.78295544 

18.36988415 

21.39271425 

23.87957538 

25.84986282 

27.31710 


28.29051532 

28.77583530 







22 

5.874266240 

10.91416887 

15.28111566 

19.04776222 

22.25706048 

24.93740208 

27.10856575 

28.7845843 


29.97529569 

30.68723324 

30.92413925 






23 

6.029268959 

11.23128639 

15.76669765 

19.70800760 

23.09824363 

25.96600770 

28.33139508 

30.20885539 


31.6087 

32.53825677 

33.00191145 






24 

6.180760221 

11.54099143 

16.24059255 

20.35190233 

23.91799356 

26.96763670 

29.52118965 

31.59345070 


33.19521851 

34.33420853 

35.01560955 

35.24241668 





25 

6.328966176 

11.84377367 

16.70359169 

20.98058126 

24.71783628 

27.94426160 

30.68043376 

32.94144468 


34.73846196 

36.07965 

36.97074938 

37.41540822 





26 

6.474089790 

12.14007146 

17.15640150 

21.59505463 

25.49912605 

28.89762693 

31.81131634 

34.25553479 


36.24176315 

37.778550 

38.87208994 

39.5265884 

39.7444875 




27 

6.616314036 

12.43027900 

17.59965570 

22.19622631 

26.26307117 

29.82928451 

32.9157779 

35.5381033 


37.70800940 

39.43437 

40.72377 

41.58088866 

42.00871 




28 

6.755804540 

12.71475230 

18.03392534 

22.78490885 

27.01075509 

30.7406219 

33.995548 

36.7912671 


39.139756 

41.05017 

42.5294331 

43.58261037 

44.21317022 

44.42314030 



29 

6.892711795 

12.99381416 

18.45972710 

23.3618358 

27.7431536 

31.63288 

35.05217 

38.016916 


40.53927926 

42.6286 

44.29227783 

45.5355 

46.36230034 

46.775 



30 

7.027173023 

13.26775837 

18.8775302 

23.9276722 

28.461149 

32.50720 

36.0870 

39.2167 


41.90861428 

44.172 

46.01515 

47.44300481 

48.46001639 

49.069 

49.271 
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